Heavy mass expansion, light-by-light scattering and the anomalous magnetic moment 

of the muon 
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Contributions from light-by-light scattering to ((? M — 2)/2, the anomalous magnetic moment of the 
muon, are mediated by the exchange of charged fermions or scalar bosons. Assuming large masses 
M for the virtual particles and employing the technique of large mass expansion, analytical results 
are obtained for virtual fermions and scalars in the form of a series in (m M /A/) 2 . This scries is 
well convergent even for the case M = m M . For virtual fermions, the expansion confirms published 
analytical formulae. For virtual scalars, the result can be used to evaluate the contribution from 
charged pions. In this case our result confirms already available numerical evaluations, however, it 
is significantly more precise. 
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High precision measurements of the muon anomalous 
magnetic moment, = (g^ — 2)/2, are used for stringent 
quantitative tests of the theories suggested for describing 
particle interactions. Presently the world average of the 
muon anomalous magnetic moment has a relative preci- 
sion of 0.7 parts per million (ppm) [l|,|2j 



4 (exp) = 11659203(8) x 10 
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and a further reduction of the experimental error by a 
factor two is within reach. 

In the standard model (SM) the theoretical value of 
a M is given by a sum of three contributions: a M (SM) = 
a^(QED) + a M (weak) + a^(had). The contribution from 
QED, a^(QED), which includes those from virtual lep- 
tons, and the one from weak interactions, a M (weak), 
can be uniquely evaluated in perturbation theory with 
the results: a M (QED) = 11658470.57(0.29) x 10" 10 
and a M (weak) = 15.1(0.4) x 10" 10 4]. The hadronic 
piece a AI (had), however, is sensitive to long distance ef- 
fects and cannot be evaluated in a perturbative frame- 
work. Using the experimental value from Eq. Q one ex- 
pects for the remaining hadronic contribution (had) = 
a^(exp) - o /1 (QED) - a M (weak) = 717(8) x 10" 10 . To 
test the standard model and to search for "new physics" , 
this value must be reproduced by a precise evaluation of 
the effects to be discussed in the following. 

The dominant contributions to a M (had) are those 
from the hadronic vacuum polarisation in lowest order 
(one particle irreducible part) with most recent results 
a^had.LO) = (702±6±14) x lO" 10 5], a M (had,LO) = 
692(6) x 10" 10 6], a M (had,LO) = 683.6(8.6) x 10" 10 0, 
a^had.LO) = (683 ±5.9 ±2.0) x 10~ 10 g and the range 
between a M (had, LO) = (684.7 ± 6.0 cxp ± 3.6 rad ) x 10~ 10 
and a M (had, LO) = (709.0 ± 5.1 cxp ± 1.2 rad ± 2.8 SU(2) ) x 
10~ 10 for analyses based on e + e~ and r data sets respec- 
tively 0. 



The next-to-leading order receives one contribution 
from the reiteration of the hadronic vacuum polarisa- 
tion which is known to be negative and can be cal- 
culated unambiguously and with sufficient precision, 
a^(had,NLO;pol) = -10.1(0.6) x 10" 10 [l<J. Its sign 
is strictly fixed: the kernel is negative definite and 
the spectral density of the two-point correlator of the 
hadronic electromagnetic current is positive. Another 
NLO contribution originates from light- by- light scatter- 
ing which cannot be calculated on the basis of data 
or first principles. The evaluations of this term have 
changed dramatically during last years 0,0,0,0,0], 
with a M (had,NLO;lbl) = 8.6(3.2) x 10" 10 as the present 
conservative estimate based on the evaluation of hadronic 
formfactors and resonance contributions 0, 0|. It is 
remarkable that the approach of ref. based on the 
concept of quark-hadron duality with an effective quark 
mass accounting for the QCD long distance effects leads 
to a M (had, NLO, lbl)d ua iity = 14(3) x 10 -10 in good agree- 
ment with the evaluations based on models for the 
hadronic formfactors. 

One important part of light-by-light-scattering ampli- 
tudes originates from neutral, low mass intermediate res- 
onances, dominantly ir° and, less important, rj, another 
from charged pion loops. This second term is the main 
subject of the present paper. 

In view of the low mass of the muon and pion, as com- 
pared to the characteristic scale of the pion form factor, 
to 2 w 0.6 GeV 2 , pions can be treated as pointlike. The 
corresponding problem of loops of (pointlike) fermions 
has been solved in analytical form, with the final result 
expressed in the form of a complicated expression of poly- 
logarithms depending on (rrij/m 2 ) 0|. 

In the present paper, a different approach is adopted, 
which makes use of the fact that m 2 /™ 2 f» 0.6 is sig- 
nificantly smaller than 1. Using the heavy mass expan- 



sion [16(, the result is obtained in the form of a power 
series, which can be systematically constructed up to ar- 
bitrary orders in (m M /m w ) 2 . The corresponding expan- 
sion can also be constructed for heavy fermions, leptons 
as well as quarks. The fermionic contribution to light- 
by-light scattering is given by the 6 diagrams depicted 
in Fig. ^ (with permutations of legs omitted). The an- 
alytical value for this contribution is known in QED for 
arbitrary values of m and M since long ^3] • Also the 




FIG. 1: QED-type light- by- light diagrams 



two leading terms in the m 2 /M 2 expansion can be found 
in this reference. Below we reproduce these terms and 
give three more terms of the expansion. To evaluate the 
contribution of charged pion loops we consider the pion as 
an elementary pointlike particle and use the Lagrangian 

i S QED = \D^\ 2 - ™ 2 M 2 , D^ = d^- ieA^ (2) 

to describe its interaction with photons. This is well jus- 
tified since the integral is convergent in the high energy 
region, and the internal structure of the pion is not yet 
resolved for momenta of order m M or m v . 

The 21 diagrams for the light-by- light contribution as 
derived from scalar QED are displayed in Fig. |3 (permu- 
tations of legs omitted) . Numerical results for this contri- 
bution were already reported in refs. 

[HEH Analytical 

expressions, however, are presented in this work for the 
first time. In view of the fact that this contribution in 
indeed comparable with the expected experimental error, 
an independent confirmation seems furthermore, highly 
desirable. 

To compute the contributions of different diagrams we 
first project onto the relevant form factor of the anoma- 
lous magnetic moment . Subsequently we employ the 
well-known method of asymptotic expansion in the small 
mass ratio m/M [l^. As an example let us consider the 
large mass expansion of the generic diagram depicted in 
Fig. [3J To obtain the asymptotic expansion for a given 
diagram one has to compute the sum of different contri- 
butions which are simpler than the original Feynman in- 
tegral. These different contributions can be classified ac- 
cording to the so called " hard" subgraphs of the original 
diagram T. These are defined as subgraphs containing all 
heavy lines such that the corresponding co-subgraph re- 
mains one-particle irreducible. The last step of the heavy 
mass expansion consists of the simple Taylor expansion 
of the "hard" subgraphs in small masses and external 
momenta. The generic prototype diagram without the 



2 




FIG. 2: Scalar QED light-by-light diagrams 

external photon line shown in Fig. [3] leads to four dif- 
ferent "hard" subgraphs: 1. the graph T itself; 2. the 
subgraph formed by the lines (1,2,3); 3. the subgraph 
formed by the lines (1,2,3,4,5,7); 4. the subgraph formed 
by the lines (1,2,3,5,6,8). We now discuss these pieces 
separately. 
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FIG. 3: A prototype diagram. Lines 1, 2 and 3 have the mass 
of pion, lines 7 and 8 — the mass of muon while lines 4, 5 and 
6 are masseless. 

1. In this case we perform a Taylor expansion of the 
diagram in Fig. Oin the external muon momentum p and 
the mass m. The resulting tensors p^p^ ■ ■ -Pfin can be 
treated with the use of the combinatoric formula 

r(d/2) , 2 

-» 2 n T(d/2 + n) {P ' 5s ' (n) 
£Vi2W--Pm2„+i -> (3) 

where g s ,(n) is the totally symmetric product of n metric 
tensors g^u, d is the space time dimension within dimen- 
sion regularization. The resulting scalar one-scale 3- loop 
bubbles where computed using the package MATAD [2fJ| . 

2. Here we perform the Taylor expansion in the mo- 
menta external to the subgraph formed by lines (1,2,3). 
The resulting symmetric tensors built from the loop mo- 
mentum of the "hard" one-loop tadpole integrals are 
expressed by the squared "hard" momentum with the 
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help of Eq. J2J • The remaining scalar one-loop tadpoles 
are expressed in terms of T- functions. Two-loop scalar 
on-shell integrals are evaluated with the package ON- 
SHELLY [m. 

3. For this subgraph and similarly for subgraph (4) we 
perform a Taylor expansion in the external momenta and 
the muon mass m. The resulting tensor one- loop on-shell 
integrals over soft momenta can be reduced to scalar in- 
tegrals for which simple explicit results exist (e.g. [13] )• 
Another way to handle the tensor structures that has 
been used for a check of the calculation is the reduction 
of the tensor integrals for a hard subgraph (explicit for- 
mulae are given in [23j|L The remaining one- loop on-shell 
and two-loop single scale tadpole integrals can be finally 
reduced to F-functions. 

The steps described above were automized with the 
Feynman diagram analyser DIANA |24| and the com- 
puter algebra system FORM [2j| . 

To discuss our results we introduce the standard nor- 
malization for the light-by-light contribution 
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For the fermionic light-by-light contribution we find 

^ 2 '3. 19 N 
16, 



aAl ( 77 ;QED) = ^ I -C 3 -- 



m 4 /13 



~Af 4 I 18^ 



161 
1620 



<2 



831931 161 



972000 3240 



-L — 



16189 
97200 1 



m 6 /17 13 1840256147 
~M6 I 36^ 3 ~~ 224^ ~~ 3556224000 



4381 
120960' 



-L 2 



-L 



24761 
317520" 
2047 



:Ca- 



453410778211 



-L 2 



5207 

189000" 
10 / 5 



54000 ~ 1200225600000 
41940853 



-L 



952560000 



m 



M io \Jg& 

328337 , 
14968800 



1187 
44550 



(2 



86251554753071 
287550049248000 

12 



640572781 \ n I m 
-L) +0 



23051952000 



M 



12 



(5) 



where L = ln(M 2 /m 2 ), m and M denoting muon and 
fermion mass respectively, C2 = C(2) = 7r 2 /6, £3 = C(3). 
The first two terms of this ex pan sion coincide with the 
result given explicitly in ref. [17]. the other terms are 
new |26j. 

The expansion in scalar QED, relevant for the charged 
pion contribution, has the form 

1 37 x 
M 2 V4 ~ 96 / 
282319 67 
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4 / 1 67 
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These results where obtained in a general covariant 
gauge, thus providing additional checks at different steps 
of the calculation. 

Let us discuss some general features of both series. The 
leading term does not contain logarithms of the mass ra- 
tio and represents therefore a pure hard contribution. It 
is obtained from the direct expansion of the graph in the 
small parameters of the problem: the muon momentum p 
and the muon mass m. The integral is infrared finite and 
no soft subgraphs appear in this order. It is represented 
by a one-scale integral with the scale given by M. In the 
next order of m 2 jM 2 soft subgraphs appear and, conse- 
quently, logarithms of the mass ratio. This is expected 
on general grounds: A gauge invariant effective action 
proportional to (F^) 4 / M 4 — the Euler-Heisenberg La- 
grangian — can be constructed which is induced by heavy 
particles. In all orders of the expansion the maximal 
power of the logarithm is two, a consequence of the sin- 
gularities of the Feynman diagrams, which is reflected in 
the structure of the expansion. 

From Eqs. (|5I6[) it is apparent that the result is sensi- 
tive to loop momenta of order M. Otherwise the leading 
order would behave as 1/M 4 and could be easily obtained 
by using the local Euler-Heisenberg Lagrangian. 

Numerically we get for the fermionic contribution (x = 
m 2 /M 2 ) 

a M ( 77 ;QED) = 0.6156x 

+ (-0.1512 + 0.1666 ln(x) - 0.04971n 2 (x))a; 2 

+ (-0.0453 + 0.0780 ln(x) - 0.0362 In 2 (x))x 3 (7) 

+ (-0.0194 + 0.0440 ln(a;) - 0.0276 ln 2 (x))a; 4 

+ (-0.0099 + 0.0278 ln(x) - 0.0219 In 2 (x))x 5 + 0(x 5 ). 

The series seems to converge even at the point x = 1 
where the sum of five terms a^{m = M) = 0.616 — 
0.151 - 0.045 - 0.019 - 0.010 . . . leads to 0.39, which has 
to be compared with the exact result a M (m = M)| oxact = 
0.3710... IHI23. 

In scalar QED the expansion has the following numer- 
ical form 

a^jr^QED) = -0.0849x 
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+ (0.0220- 0.0194 ln(ar) +0.0052 ln 2 (x))a; 2 

+ (0.0049 - 0.0067 ln(x) + 0.0027 In 2 (x))x 3 (8) 

+ (0.0016 - 0.0029 ln(x) + 0.0015 ln 2 (x))x 4 

+ (0.0007 - 0.0015 ln(x) +0.0009 ln 2 (x))x 5 +0{x 5 ). 

At the point m = M one finds 0^(77; sQED) = 
-0.0849 + 0.0220 + 0.0049 + 0.0016 + 0.0007 = -0.0557 
with reasonably decreasing individual terms. In this case 
the exact result is unknown, it could, however be cal- 
culated with the technique of ref. [17j. The difference 
between the contributions of fcrmions and scalars is ap- 
parent: for scalars the result is negative and almost an 
order of magnitude smaller in absolute value. The rate 
of convergence is quite similar in both cases. 

In the present case of interest the actual value of the 
expansion parameter is given by the ratio (m^/m^) 2 = 
(106/139. 6) 2 = 0.577. This leads to a rapidly converging 
series 0^(77; sQED) = -0.0490(1-0.233-0.037-0.008- 
0.002) = —0.0353 with an estimated accuracy better than 
0.5%. This value should be compared with the numeri- 
cal evaluation of ref. where two methods were used 
with the results 0^(77; sQED)| numi i = —0.0437(36) and 
a M (77;sQED)| num , 2 = -0.0383(20). The slight disagree- 
ment between these two values is presumably caused by, 
quoting the authors of ref. .18], "a delicateness of can- 
cellation when separate contributions are put together" . 
The present approach is based on an analytical evalua- 
tion that guarantees the absence of error accumulation. 
The accuracy of the result is determined by the conver- 
gence rate of the series in Eq. © which is quite high for 
physical values of muon and pion masses. 

To summarize: A systematic expansion of the light-by- 
light contribution to the anomalous magnetic moment of 
the muon has been presented which is valid for virtual 
fermions and bosons with masses above or equal m^. 
The formulae are sufficiently accurate for the physical 
applications of interest. 
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